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LHC Implications of the WIMP Miracle and Grand Unification 
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With the assumptions that dark matter consists of an electroweak triplet and tliat the gauge 
couplings unify at a high scale, we identify robust phenomenological trends of possible matter 
contents at the TeV scale. In particular, we expect new colored states within the LHC reach that 
can have Yukawa couplings A to quarks and the Higgs. We investigate the collider signatures that 
are characteristic of all such models by adopting the model with the simplest matter content as 
a benchmark. The A couplings are constrained by flavor/CP physics. In the largest portion of 
the allowed parameter space the new colored particles are stable on collider time scales, hence 
appearing as R-hadrons, for which there is discovery potential at the early LHC {^/s = 7 TeV, 
1 fb~^). Flavor/CP constraints nevertheless do allow a sizable range of A where the new colored 
particles decay promptly, providing a new Higgs production channel with a cross-section governed 
by the strong interaction. Studying the case of ft —> WW, we show that it is possible for the Higgs 
production from this new channel to be discovered before that from the Standard Model at the 
LHC. 



I. INTRODUCTION 

The existence of dark matter (DM) is arguably the 
most compelling evidence for new physics beyond the 
standard model (SM). Even though existing data pro- 
vides little insight into the identity and nature of the 
DM particle, a simple and robust candidate is provided 
by a weakly interacting massive particle (WIMP), as its 
relic abundance will be automatically of the right size 
when its mass is at the TcV scale ("the WIMP miracle"). 
Since the null results from direct DM search experiments 
have excluded a WIMP with nonzero hypercharge [l| , the 
simplest WIMP DM candidate is an SU(2)L-triplet with 

Y = 0, denoted as a V hereafter. The V can be made 
stable by imposing e.g. the Z2 parity V —V. A care- 
ful calculation of the relic abundance of a 1^ is performed 
in Ref. [2[, including non-perturbative effects and possi- 
ble co-annihilations, finding that the V mass should be 
2.5 TeV if the V is spin-0, or 2.7 TeV if spin-1/2, assum- 
ing that the V accounts for the entire missing mass of 
the universe. 

Unfortunately, since the V is heavy and color-neutral, 
it is virtually impossible to be directly produced from 
collisions at the Large Hadron Collider (LHC). However, 

V may well be part of a bigger, well-motivated extension 
of the SM containing other new particles that can give 
rise to observable LHC signals. 

In this paper, we adopt gauge coupling unification [sj 
as a guiding principle besides WIMP DM. Although this 
is still not constraining enough to point to one single 
model, it is possible to identify generic, robust phe- 
nomenological trends in such extensions of the SM. Con- 
cretely, we adopt the as a dark matter candidate and 
demand perturbative gauge coupling unification.^ We 



assume no extra mass scales other than the unification 
scale and the TeV scale dictated by the WIMP mira- 
cle. Finally, we assume that new particles, including the 
V , are all fermions to avoid extra fine-tuning problems 
associated with scalar masses. These assumptions read- 
ily imply the existence of additional new particles at the 
TeV scale, since the gauge couplings in the SM-I-F the- 
ory do not unify. In particular, we find that there must 
exist new colored particles at the TeV scale. We also find 
that the new colored particles generically allow Yukawa 
couplings to quarks and the Higgs. 

Split supersymmetry [1, Q is a well-studied scenario 
also based on WIMP DM and gauge coupling unifica- 
tion, as well as the absence of supersymmetry at the TeV 
scale as in our scenario. There are two major differences. 
First, contrary to one of our assumptions above, split 
supersymmetry has an extra threshold between the TeV 
and unification scales, where we have all the squarks and 
sleptons, and most importantly, the second Higgs dou- 
blet, which affects unification. Second, in split supersym- 
metry, WIMP DM has to be a nontrivial composition of 
the higgsinos, wino and/or bino; if we assume that dark 
matter in split supersymmetry is a pure wino (i.e. the 
V) as in our scenario, gauge coupling unification would 
not work well (using the same criteria for precision as 
we use in Sec. |ll| unless we have a hierarchy larger than 
two orders-of-magnitude between the higgsino and gluino 
masses. 

The existence of new colored particles with Yukawa 
couplings A to quarks and the Higgs suggests the fol- 
lowing two scenarios for the LHC. If A is sufficiently 
small, the new colored particles will be collider stable, ap- 
pearing as massive stable hadrons ("R-hadrons"). Since 
R-hadron signals can be quite spectacular, this is an 
exciting possibility already for the early LHC run at 



^ Note that a y is not only the simplest WIMP DM candidate 
of all but also the only SU(2)l multiplet with zero hypercharge 



that appears within simple SU(5) multiplcts, 5, 10, 15, and 24. 
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^/s = 7 TcV with an integrated luminosity 1 fb""'^. 
On the other hand, if A is not so small, the new col- 
ored particles will decay promptly via A, with an 0(1) 
fraction of their decays containing Higgs bosons. This is 
an interesting new production channel for the Higgs bo- 
son, where the size of the cross-section is governed by the 
strong interaction, potentially making the LHC a "Higgs 
factory" . 

To perform quantitative benchmark studies of these 
characteristic phenomenological features of WIMP DM 
and unification, we choose as a simple benchmark model 
consisting of a DM candidate V and new colored parti- 
cles X (to be specified more explicitly later). Among all 
models with WIMP DM and unification, this benchmark 
model contains the smallest number of new multiplets 
beyond the SM, but it already exhibits the two classes of 
generic collider signatures mentioned above. 

Our analysis on this benchmark model will show that, 
for the range mx = 360-650 GeV, the early LHC phase 
(7 TcV; 1 fb^^) should have sufficient discovery potential 
for the R-hadron case. For the Higgs factory case, we 
will lay out an experimental strategy for the full LHC at 
^/s = 14 TeV. This consists of two parts, the discovery 
of the X and measurement of mx, and the discovery of 
the Higgs bosons from the X decays. We will show that 
with 10 fb"^ of data at the LHC (14 TeV), it should be 
possible to discover the X and the Higgs bosons from 
the X decays in the range 300 < mx ^ 550 GeV for 
a moderately heavy Higgs (i.e. decaying to weak gauge 
bosons). 

This paper is organized as follows. In Sec.|ll]we survey 
possible extensions of the SM that feature gauge coupling 
unification and WIMP DM. In Sec. IIIII we describe our 
simple benchmark model that phenomenologically repre- 
sents all such extensions. The couplings A of the new 
colored particles to quarks and the Higgs are expected to 
have an upper bound from flavor/CP/electroweak con- 
straints, which is analyzed in Sec. IIVI using the bench- 
mark model. The collider signatures of the R-Hadron 
case and the Higgs factory case are studied in detail in 
Sec. IV Al and Sec. IV Bl respectively. In Sec. |Vl]we sum- 
marize our analyses and indicate some possible future 
directions. In Appendix |X] we discuss how proton decay 
can be avoided in the class of models we consider in this 
paper. In Appendix IB] we comment that the addition of 
higher-dimensional operators to our Lagrangian does not 
have any impact on our analysis. 



II. EXTENSIONS OF THE SM FEATURING 
WIMP DM AND COUPLING UNIFICATION 

In this section we enumerate possible extensions of the 
SM that contain the WIMP DM candidate V and are 
consistent with gauge coupling unification, and identify 
generic features shared by such extensions. This analysis 
will serve as a basis for our choice of a benchmark model 
in Sec. imi 



Let us assume unification of the SM gauge group into a 
simple group (such as SU(5)) to fix the normalization of 
hypercharge. The DM candidate V can be embedded into 
a 24 of SU(5), consistent with this assumption. More- 
over, let us assume that all new particles, including the 
V, are spin- 1/2 fermions in order to avoid unnecessary 
extra fine-tuning problems associated with scalar masses 
besides the notorious existing problem with the Higgs 
mass. As calculated in Ref. the fermionic V mass is 
fixed by the relic abundance to be mv = 2.7 TeV, which 
we will assume to be the case hereafter. 

This cannot be the end of the story, however, because 
the SM augmented by only the V is not consistent with 
gauge coupling unification. There must be additional 
new particles. In principle, these additional particles 
could appear anywhere below the unification scale. How- 
ever, since we must presume some underlying dynamics 
that generates the TeV scale in order for the WIMP mir- 
acle not to be a mere coincidence, we adopt the sim- 
plest assumption that the same dynamics also provides 
TeV-scale masses to the additional new particles, with 
no extra mass scale other than the TeV scale and the 
unification scale. We take the new fermions to be vector- 
like, because chiral fermions would require electroweak 
symmetry breaking to acquire TeV-scale masses, which 
would generically lead to dangerously large corrections 
to precision electroweak observables, in particular the p 
parameter Q. Let us further restrict ourselves to the 
case where the vectorlike fermions can be embedded into 
the simplest SU(5) multiplets, 5 © 5, 10 © TO, 15 © T5 
and 24. Thus, we consider 

Q ^ (3, 2)1/6 , U-^{^, 1)2/3 , ^ ~ (3, l)_i/3 , 
L ^ (1, 2)_i/2 , S ^ (1, l)_i , A ~ (3, 2)_5/6 , 
5^ (6,l)_2/3, T~ (1,3)1, (1,3)0, 
G~(8,l)o, 

as well as the conjugates Q'^, W^, • • • , T"^, except for V 
and G, which are real. (Our convention is such that H 
has the quantum numbers (1, 2)1/2. The SM fermions 
are denoted by lower-case letters, q, u^, (f , £, e^.) As we 
will see below, this already provides a sufficient number 
of candidate models for us to observe generic trends in 
extensions of the SM with WIMP DM and gauge coupling 
unification. 

In searching for possible field contents that can lead to 
unification, there are various uncertainties that must be 
taken into account when we "predict" the SU(3)c cou- 
pling as in terms of ai_2! which we regard as precise. 
First, in our RG analysis, which we perform at the 1-loop 
level, there is a threshold ambiguity at my = 2.7 TcV. 
We estimate this uncertainty by varying the MS subtrac- 
tion scale /i from my /V2 to \/2 my- Second, unlike the 
V mass, the masses of the additional fermions cannot 
be fixed a priori and can be anywhere at the TeV scale, 
from a few hundred GeV to several TeV. Therefore, we 
scan over the additional fermions' masses in the range 
between 300 GeV and 10 TeV, where for simplicity we 
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2X 2X + D + L 2U + T + E 

U + S + 2T 2X + 2f/ + T3X + L + £ 
3X + G + y 3C/ + T + -L 



TABLE I: Possible combinations of new fermions that, to- 
gether with the DM V , could lead to unification. It is under- 
stood that the new fermions are vectorlike, so in writing X, 
[/,••■, the presence of their charge conjugates X'^ , , ■ ■ ■ is 
also implied, except for the Majorana fermions V and G. 

assume a single common mass for all of them.^ These 
two sources of uncertainty each shift our prediction 
at the level of a few times iS.a"^^'^^ , the experimental un- 
certainty in the measurement a^^''^^ = 0.1184±0.0007 0]. 
We demand that the "band" in our as prediction com- 
bining these two uncertainties have an overlap with the 
band corresponding to BAckg'^^^'' (i.e. Scr). There are also 
threshold effects from unspecified GUT physics, but we 
assume that they are similar in size to the uncertainties 
mentioned above and simply neglect them. Finally, we 
demand that the coupling at the unification scale be per- 
turbative, acuT < 1- 

A similar analysis was performed in Ref. where 
the main difference lies in the treatment of proton decay. 
While Ref. Q demands the unification scale to be higher 
than ~ 10^^ GeV in order to sufficiently suppress proton 
decay, we choose to impose a symmetry to forbid proton 
decay, and only demand the GUT scale to be higher than 
10^ TeV (and lower than the Planck scale ~ 10^® GeV) to 
avoid having to address possible confiicts between GUT 
physics and flavor/CP bounds. The use of a symmetry 
to forbid proton decay requires some model building at 
the unification scale, but it has no observational conse- 
quences for the TeV-scale physics, so we leave the model 
building to Appendix [X] Another difference between our 
analysis and that of Ref. @ is that in calculating the run- 
ning of the gauge couplings, Ref. Q assumes that all new 
particles have masses near mz^ while our masses span a 
wide range around a TeV, as we described above. 

There are 22 models satisfying the above criteria with 
no more than 3 types of multiplets in addition to the V 
and no more than 3 generations per type. Particularly 
simple ones are listed in Table HI All the 22 models share 
the following features: 

(1) There exist new colored particles. 

(2) The quantum numbers of these new colored par- 
ticles allow Yukawa couplings to quarks and the 
Higgs. 

Property (1) is clearly favorable for hadron colliders. 
Even better, all the 22 models actually survive even if 



^ New charged/colored particles below 300 GeV arc likely to be 
already excluded, as will be illustrated by the analysis of the 
benchmark model below. 



we restrict the additional fermions' masses to the range 
300 GeV-1 TeV, so they all can be potentially within 
the LHC reach. Property (2) is not satisfied by S and 
G, but they only appear twice each among the 22 mod- 
els, and even those models contain other colored particles 
that do satisfy property (2). We therefore identify these 
properties as robust LHC implications of WIMP DM and 
unification.'^ 

To assess the robustness of the above features further, 
one can repeat the exercise with more conservative esti- 
mates on the uncertainties in the prediction of a^. For 
example, if we vary the matching scale /i from my/2 to 
2'mv (with everything else treated as above), we obtain 
63 models, of which there is only one model {V-\-E) with- 
out colored particles,^ and only 6 colored models without 
property (2). Again, most models survive even if we re- 
strict the search in the "LHC-accessible" range 300 GeV- 
1 TeV; 44 models in total, only one colorless model, and 
only one colorful model without property (2). 

In the next section, we choose a benchmark model that 
represents characteristic phenomenologies of all these 
models which follow from properties (1) and (2). We will 
then use the benchmark model for further, more quanti- 
tative analyses in later sections. 



III. THE BENCHMARK MODEL 

Given the insights from Sec. |IT1 we select the model 
with the dark matter V and two generations of X © X"^ 
as our benchmark. This is the simplest of all models with 
a V featuring unification, having the smallest number of 
new multiplets beyond the SM. But most importantly, 
the collider phenomenology of this model is representa- 
tive of all models identified in Sec. HI] as far as the LHC 
phenomenology is concerned. 

The most general renormalizable Lagrangian for the 
2X + V model consistent with the Z2 symmetry V ^ ~V 
reads 

+ J2 [^afT^iD^Aa + A^^ID^A:^ (1) 

a 

-{inx^X^X^ + J2 A- Hd'rX^ + C.C.J , 

i 

where a = 1, 2 and i = 1,2,3 denote generations of the 
A and d-type quark, respectively. We also use A_i/3 



^ Recall the crucial role of the WIMP miracle in selecting the TeV 

scale as the mass scale for the new particles. 
^ Actually, a closer inspection reveals that unification favors the E 

in the V + E model to be lighter than fa 300 GeV. The quantum 

numbers of the E allow it to decay to Z in particular, so this 

model is already excluded. 
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and ^-4/3 to refer to the upper and lower SU(2)l com- 
ponents of X, respectively, where the subscripts denote 
the electric charges. 

It is technically natural for A to take any value, but 
there are obvious phenomenological constraints. First, A 
has to be nonzero because the Xs must eventually de- 
cay to avoid cosmological problems. (However, the Xs 
could decay via higher dimensional operators. See Ap- 
pendix [B]). Second, A must be much less than 0(1) be- 
cause A breaks the U(3)^ flavor symmetry of the SM, pro- 
viding new sources of flavor/CP violations in addition to 
the SM Yukawa couplings. We will return to flavor/CP 
constraints in Sec. IIVI 

The leading decays of X induced by A can be most 
easily understood by the Goldstone equivalence theorem. 
In the limit of keeping only the X mass, the equivalence 
theorem tells us that X_i/^ will decay as 

X_i/3^Z + d, X_i/3^/l + d (2) 

with equal probabilities, where d can be any down-type 
quark. ^ Note that the equivalence theorem holds only 
in this limit. When the finite masses of the Higgs and 
Z bosons are taken into account, the branching fraction 
for the Z channel is expected to be somewhat larger due 
to phase space. For low X masses and/or a heavy Higgs, 
this can have a significant impact on the phenomenology. 
Similarly, X_4/3 can decay as 

X_4/3-> W^~+rf, (3) 

which will be the dominant decay mode as long as the cor- 
responding rate can be regarded as prompt on the collider 
time scale. When this rate drops below the displaced- 
vertex range, the dominant decay of the X_4/3 will be 
through the weak interaction to an X_i/3, which be- 
comes slightly lighter than ^-4/3 after electroweak sym- 
metry breaking, as we will elaborate more in Sec. IV Al 

Therefore, depending on the size of A, we have either 
of the following collider signatures: 

(A) If A is sufficiently tiny, the X will be stable 
on collider time scales, and upon production it 
will hadronize into stable massive hadrons ("R- 
hadrons"). R- hadrons are easy to observe when 
they are charged, so X may be discoverable already 
in the early LHC run (i.e. 7 TeV, 1 fb"^). 

(B) If A is not so small (but small enough to satisfy 
flavor/CP constraints), the X will decay within 
the detector, and as we have seen above, roughly 
a quarter of Xs (a half of X_i/^ s) will decay to a 



^ The small violation of the equivalence theorem induces additional 
decays such as X_i/^ — >■ W~ + u, which can be thought of as 
arising from mixing of the Xs with down-type quarks. How- 
ever, as we will see in Sec. IIVI flavor/CP bounds constrain such 
mixings to be tiny, rendering these decay modes negligible. 



Higgs boson (plus a jet). This is an exciting pos- 
sibility — a "Higgs factory" — where the Higgs 
bosons are produced with a characteristic 2 — > 2 
cross-section of the strong interaction. In the re- 
maining 3/4 of the time, the X will decay to a Z 
or W boson. Then, leptonic Z decays can be used 
to discover the X itself. 

In Sec. lIVI we will show that flavor/CP constraints indeed 
allow a window for the case (B). 

Note that possibilities (A) and (B) are common to all 
the 22 models identified in Sec.|TTl For example, in mod- 
els containing Q instead of X, the A coupling in Eq. ^ 
should be replaced by XuHu'^Q + XdH*cPQ, which would 
exhibit the same phenomenology as above. In models 
containing U or D instead of X, the A coupling is re- 
placed by XHW^q and XH*D'^q, respectively, which is 
again phenomenologically equivalent. 

Actually, the models with Q, U and/or D have ad- 
ditional potentially interesting modes Q — > Z -|- i, Q — > 
W+t, U Z+t, or D ^ W+t. While the appearance of 
the top constitutes a qualitative difference in the collider 
phenomenology, a full analysis for such decay channels is 
more complicated due to the higher final-state multiplic- 
ity from the top decay. There are reducible background 
sources which cannot be simulated reliably at the matrix 
element level due to the large number of final-state par- 
ticles, and even for irreducible backgrounds the issue of 
combinatoric backgrounds makes searches more difficult. 
Some of these problems may be ameliorated if an 0(1) 
fraction of tops is produced with large pt, such that the 
recently developed methods of boosted top-tagging Q 
can be applied. We will leave these more complicated 
cases to future work and focus in this paper on the phe- 
nomenology of signatures (A) and (B). 

IV. FLAVOR/CP AND ELECTROWEAK 
CONSTRAINTS 

In this section, we analyze flavor/CP violations as well 
as corrections to precision electroweak observables in the 
benchmark model. These corrections arise due to the 
coupling 

C D A,a Hd^^Xa + c.c. , (4) 

where a = 1,2 and i ~ 1,2,3. We will adopt the most 
conservative assumption that A is an "anarchic" matrix 
without any special texture or alignment: 



A ^ 


- A 


A ^ 


- A 


A ' 


- A 



Therefore, the bounds discussed in this section could 
be relaxed by further model building or extra assump- 
tions on the structure of A. (For example, Ref. ^ in- 
troduces a model with a single X particle (and no V) 
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FIG. 1: Diagram leading to the operator ((Gj. 
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FIG. 2: Diagram leading to the four fermion operator 



with a similar coupling selectively to the third generation, 
and consequently their model is much less constrained by 
flavor/CP.) 

Since the Xs are heavier than all SM particles and 
wc anticipate a small A, we integrate out the Xs and 
analyze effective operators in powers of X/mx- Strictly 
speaking, the ratio {H) /mx is not a very small number, 
so contributions higher order in {H) /mx can change our 
estimates by an 0(1) factor. However, our interest is 
to show that a robust "Higgs factory" window can exist 
for the broad scenario of extending the SM with WIMP 
DM and unification, rather than placing precise bounds 
on this particular benchmark model. Therefore, order- 
of-magnitude estimates suffice for this purpose. 

The most stringent bound comes from K^-K^ mix- 
ing. The relevant tree and 1-loop diagrams are shown in 
Figs. [T] and [U respectively. Upon integrating out X in 
Fig. [TJ we generate the operator 



E 



A,-„A 



la a J 

57 



(6) 



Below the Z mass, using this operator twice with all the 
four Hs put to VEVs would generate 4- fermion operators 
with four right-handed down-type quarks with a coef- 
ficient ^ X^v'^ /rrix, which is conservatively ~ A^/m^. 
In particular, the imaginary part of the coefficient of 
{cP(jf^s'^){(F(j^s^ is constrained and has to be less than 
(10^ TeV)~^ With our assumption of anarchic A, 

we expect an 0(1) phase in A^, so we obtain the bound 



A < 10" 



mx 
1 TeV 



(7) 



Similarly, upon integrating out X in Fig. [21 we generate 
the operator 

^ 1^ {d^(Tf,d^)[qk(T' qi). (8) 



a,i,j,k,£ 



Xa di 



FIG. 3: One of the diagrams leading to the dipole operator 



From this, the imaginary part of the coefficient of 
{(F a ^s^){da^^ s) can be estimated to be 



1 A^flS 



IGtt^ m\ 



12 



(9) 



el w el with the Cab- 



where wc have used {y\jj^) 
bibo angle 0c ~ 1/5. This coefficient should be less than 
(10^ TeV)-2 [3, which implies 



A < 10" 



1 TeV 



(10) 



The bound PU)) is slightly stronger than the bound ([7]) 
for mx < 1 TeV, so we adopt Eq. pO)) as our upper 
bound on A. 

No other constrains are as strong as Eq. (fTO|). For 
example, let us look at the dipole operators generated 
from diagrams as in Fig. [3j 



9F 

167r2 ^ 

a,i.,j,k 



(11) 



where F ~ B,W denotes the electroweak gauge fields, 
which in particular contribute to 5 — )■ S7 after elec- 
troweak symmetry breaking. Since these operators are 
suppressed by the small bottom Yukawa coupling yi, ~ 
1/40, one can think of them as the b s-f dipole op- 
erator in minimal flavor violation [Tl] | with the scale 
A ^ Anecmx /{V^ X), where e is the QED coupling. 
Then, the bound A ~ 10 TeV from 6 ^ 57 [ul Hi im- 
plies X/mx < 10-iTeV"\ which is again weaker than 
Eq. (Uni). 

Let us also look at precision electroweak constraints. 
First, the operator ([6]) modifies Z — >■ bb: 



Zbb 



A2 



9zbb 



m 



(12) 



X 



Again, given Eq. (|10p . this is safely below the experi- 
mental bound. ^ Second, Fig. 2] generates an operator 



® Ref. 01 discusses a model with a single generation of X (and 
without V), in which they exploit this shift in Z bb to improve 
precision electroweak fits of the SM. 
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FIG. 4: Diagram contributing to the p parameter (jlSp . 



contributing the p parameter 



(13) 



The coefficient should be less then ^ 10 /v^ [Ijl with 
V ~ 174 GcV, implying A < v^mx/(l TeV), which again 
is much weaker than Eq. ([TU)) . 

Given the bound (fTO)) . there is clearly a robust "Higgs 
factory" window where the Xs decay promptly. Keeping 
only the X mass for simplicity, the decay rate of an X 
via coupling (|4]) is given by Fx ~ A^mx/167r. Therefore, 
demanding F^^^ < 10"^^ s (i.e. the decay length shorter 
than 0(0.1) mm), wc obtain the Higgs-factory window 



10" 



1 TcV 



mx 



< A < 10" 



1 TeV 



(14) 



Below the Higgs-factory window, there is a two or 
three orders-of-magnitude window where an X decays 
with a displaced vertex or in the LHC detector but with 
a macroscopic decay length. At a hadron collider, events 
where new physics only manifests itself at a macroscopic 
distance from the interaction point but still within the 
volume of the detector are challenging. Such signatures 
also appear in supersymmetric theories [l3, [I3I, hid- 
den valley models [l^, quirk models [13, and vector- 
like confinement models [l^. As discussed in section 
16 of Ref. [1^ (and references therein), when particles 
with macroscopic decay lengths arc produced, trigger ef- 
ficiencies can become a concern. Recently, the prospects 
of a similar final state from late-decaying neutralinos in 
a supersymmetric model have been carefully studied in 
Ref. [l^l, with encouraging results. While the case of 
macroscopic decay lengths is interesting, the subtleties 
involved require a study with a more sophisticated detec- 
tor simulation than PGS [2l| (which we use in our analy- 
sis). In addition, it is possible for R-hadrons to stop in 
the detector [2^ . These are important questions worth 
investigating in detail, which wc leave to future work. 



V. COLLIDER PHENOMENOLOGY 

In this section we will investigate in detail the two 
characteristic collider signatures of our scenario using 
the benchmark model. When the Xs arc collider sta- 
ble, we will show that the discovery of X is possible up 
to mx = 650 GeV in the early 7-TeV run with 1 fb"^ 
of data, and well past 1 TeV with the 14- TeV running. 
When the Xs decay promptly, we will concentrate on the 
discovery of the X and the Higgs bosons from X decays 
for early 14- TeV running (10 fb^^). In order to keep the 
analysis simple, we will restrict ourselves to the case of a 
moderately heavy Higgs {rrih = 200 GeV) , but we expect 
the Higgs discovery potential to be similarly enhanced 
for a light Higgs as well. 

Even though the benchmark model contains two gen- 
erations of X, we have no reason to expect that they 
would be exactly degenerate in mass. Since the produc- 
tion cross-section will be dominated by the lighter X, 
we will base our collider analysis on one generation of X 
only. This is the most conservative choice; if the Xs do 
happen to be nearly degenerate, this would significantly 
enhance the results below. 



A. R-hadron Signals at the LHC 

When the coupling A is sufficiently small, the X does 
not decay within the detector. The signature therefore 
is that of an "R-hadron" , that is, the X hadronizes with 
light colored degrees of freedom and the color-neutral 
bound state behaves as a (possibly charged) stable mas- 
sive particle. Before we present a quantitative analysis, 
let us dwell on a few qualitative features of this signature. 

4/3 



Firstly, note that even when A is very small, X_ 



still decays to X_i/^ via the weak interactions, and in 
terms of coUider signatures, the production of ^-4/3 is 
indistinguishable from that of X_i/3 because the decay 
products are virtually unobservable. As worked out in 
detail in Ref. [2^, X_4^/^ is expected to be heavier than 
A'_i/3 by only Amx — 0.60 GcV. The dominant decay 



mode is X_ 



4/3 



X_ 



with a partial width 



1/3 



TT through an ofF-shell W 



F =1.3 mm 



0.60 GeVV 
Amx J 



\ 



1 - 



(0.60 GcV)^ 



1 - 



(15) 



Am'i 



The smallness of the mass gap makes the tt very soft 
and thus unobservable at the LHC, and other subdom- 
inant decay modes have the same problem. Therefore, 
in analyzing the discovery potential or checking existing 
bounds, the production cross-section of Ar_4/3 should be 
added to that of X^i/-^. 

The charge of an R-hadron is crucial for prospects of 
observing it. In particular, the most effective way to 
trigger on a stable charged massive particle is via the 
muon system [g^]- The bound states of an X_i/3 can 
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be mesons (Xq) or baryons (Xqq). The physics of these 
bound states can be understood by regarding the X_i/3 
as a heavy version of the b quark, which is already much 
heavier than Aqcd- The hghtest B mesons are and 
B^ , with only a few-hundred-keV mass splitting, while 
the lightest B baryon A° is heavier than the S"'^ by 
340 MeV. Therefore, we expect that the lightest X- 
mcson should be lighter than the lightest X-baryon also 
by ^ 340 MeV, with a few- hundred- keV mass splitting 
between the neutral and charged X-mesons. Since the 
splitting between the lightest X-meson and X-baryon is 
on the order of Aqcd itself, we expect that an X should 
preferentially hadronize into an X-meson, which can be 
either charged or neutral with 50% probability because 
their mass difference is tinyJ 

In order to estimate the trigger efficiencies, we will use 
the following assumptions in the rest of our analysis: 

• There is a 50% chance that an R-hadron is charged 
when produced at the primary interaction. 

• This charge is retained until the calorimeter is 
reached. 

• One or more charge exchange interactions take 
place in the calorimeter, randomizing the charge of 
the R-hadron such that there is a 50% chance that 
it reaches the muon chamber as a charged particle.^ 

One of the requirements for triggering is that the particle 
reaches the muon system with nonzero charge. Most ex- 
perimental searches for massive stable particles also use 
as a selection criterion that there should be a charged 
track in the inner part of the detector that matches the 
hit in the muon chamber, even if this is not required for 



The few-hundred-keV mass difference between the {X_i/^d) me- 
son and the {X_i/^u) meson might allow one to decay to the 
other via the weak interaction, if the mass difference is larger 
than the electron mass. But such a decay would occur with an 
extremely long lifetime (like the /3-decay of the neutron), so it 
can be ignored on collider time scales. 
* In addition to conversion between the charged and neutral X- 
mesons, for which there is a tiny energetic cost of a few hundred 
keV, there is also a process where an Jf-meson scatters into an X- 
baryon in the calorimeter [25l |. which, however, requires an energy 
of at least ~ 340 MeV. Using the analogy with the b system, the 
lightest X-baryon (analogous to AjJ) should be neutral, while the 
lightest charged X-baryon (analogous to S^) should be heavier 
by ~ 190 MeV. So, a charged X-baryon, even if produced, would 
promptly decay to a neutral X-baryon (by emitting a pion) which 
would then not be caught by the muon chamber, potentially 
hurting our R-hadron signal. The question of how frequently 
this meson-to-baryon conversion occurs is highly nontrivial and 
beyond the scope of this paper. However, note that a meson- 
to-baryon conversion would not occur to the X due to the lack 
of anti-nucleons in the calorimeter. Therefore, even in the worst 
case where we always lose the R-hadrons from the X , we still have 
those from the X , so the effective cross-section would be roughly 
halved, which would correspond a small shift {r^ 0(10) GeV) in 
the mass scale. 
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FIG. 5: The cross-section for R-hadron production at the 
Tevatron after all selection cuts. The bound from the 
CHAMP search is included for comparison. 



triggering. Therefore we will also adopt this as one of 
our event selection criteria. 

Finally, an important kinematic variable, especially at 
the LHC (where the detectors are physically larger and 
the time between bunch crossings is short), is the "time- 
lag". This is defined as how much later the massive R- 
hadron (/3 < 1) arrives at the muon chamber compared 
to a relativistic particle (e.g. a muon). To be conserva- 
tive, we use the physical dimensions of the ATLAS detec- 
tor (which is larger than CMS) as described in Ref. j23 |. 
Specifically, we differentiate the barrel region {\r]\ < 1.4) 
from the end-cap region (1.4 < \r]\ < 2.5). For the bar- 
rel region we calculate the time to get to a radius of 
7.5 meters from the interaction point, and for the end- 
cap we calculate the time to get to |z| = 14.5 meters. For 
the LHC search, we will use as one of the event selection 
criteria that at least one of the R-hadrons reaches the 
muon chamber before relativistic particles from the next 
bunch crossing, i.e., with a time-lag of less than 25 ns. 

Before we investigate the discovery potential of X at 
the early LHC, let us address constraints on R-hadron 
production from the Tevatron [2^ . We will use the 
event selection criteria described in Ref. in order to 
estimate acceptance and trigger rates in the benchmark 
model. In particular we demand that events contain at 
least one particle that has \t]\ < 0.7, pt > 40 GeV and 
0.4 < f3 < 0.9, leaves a track, and is charged when it 
arrives at the muon system. For a single R-hadron sat- 
isfying these cuts, our above assumptions on the charges 
of R-hadrons imply a 25% probability for being detected. 
However, when both of the pair-produced R-hadrons are 
within acceptance and charged throughout the detector 
(a 1/16 probability), we need to correct for the fact that 
the reconstruction and trigger efficiency used in the anal- 
ysis of Ref. [1^ applies to a single R-hadron. Therefore 
the "effective cross section" for such events needs to be 
multiplied by a correction factor ^ before the compari- 
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FIG. 6: The rapidity distribution of the R-hadrons at the 
LHC (7 TeV) for the three mass points. 



FIG. 7: Time-lag distribution of the first R-hadron to arrive 
in the muon system at the LHC (7 TeV). 



son with the bound of Ref. [26|. This correction factor 
is given by ^ = (1 — (1 — e)^)/e, where e = 0.533 is the 
reconstruction efficiency for a single R-hadron We 
use CalcHEP 2 . 5 . 4 [111 with CTEQ6 parton distribution 
functions [s^l to simulate the parton level process. The 
effective cross-section after the selection cuts is plotted in 
Fig. [S] against the bound of Ref. [1^. We conclude that 
mx > 360 GeV is not excluded. 

We now turn to production of R-hadrons at the early 
LHC, i.e., for the 7- TeV running with 1 fb~^ of inte- 
grated luminosity. We choose to look at mass points 
mx = 450 GeV, 550 GeV and 650 GeV. In Fig. M we 
plot the rapidity distribution of the R-hadrons for each 
mass point to show that R-hadron production is dom- 
inantly central. As event selection criteria, we impose 
that at least one R-hadron must reach the muon detector 
(1 77 1 < 2.5) with pt > 30 GeV and a time lag of less than 
25 ns. We further demand that the R-hadron in question 
leaves a track and is charged when it gets to the muon 
chamber (a 25% probability per R-hadron as before) . In 
events where both R-hadrons reach the muon chamber 
with pt > 30 GeV, we plot the time-lag of the earlier 
(later) R-hadron in Fig. [7] (Fig. [5]). Folding in the time- 
lag cut and the probability of being charged (but leav- 
ing out reconstruction efficiencies, which are unknown 
at this time), we plot the effective cross-section after all 
selection cuts in Fig. [S) Note that the efficiency of the 
selection cuts has a slightly decreasing trend at higher 
mass, because the production occurs closer to threshold 
and fewer events satisfy the time-lag cut. Requiring at 
least 10 events for discovery, we see that X masses of up 
to 650 GeV should be within reach with 1 fb~^ of data 
from the 7- TeV run. 

Using the same selection criteria for the 14- TeV run- 
ning, we also plot in Fig. [10] the effective cross-section 
after cuts in that case. We see that even with early 14- 
TeV data (e.g. ss 10 fb"^), X masses weU past 1 TeV 
should be within reach. 
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FIG. 8: Time-lag distribution of the second R-hadron to ar- 
rive in the muon system at the LHC (7 TeV). 
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FIG. 9: Total production cross-section of R-hadron produc- 
tion and effective cross-section after selection cuts at the LHC 
(7 TeV) for the three mass points. 
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FIG. 10: Effective R-hadron production cross-section after 
selection cuts at the LHC (14 TeV). 



B. The LHC as a Higgs factory 



For A in the range given by Eq. ([T4| . the X decays 
promptly. For the X_i/2^ the dominant decay modes 
are Z + j and h + j. Since production proceeds through 
QCD, this gives rise to a "Higgs factory" if mx is not too 
large. The X_4/^ decays to W + j and therefore gives no 
additional contribution to Higgs production. For the rest 
of this section, we will be interested in the production of 
X_i/3 only. We will show that there is a window in mx, 
where with 10 fb~^ of 14- TeV running, both the X and 
the Higgs can be discovered. 

We will first dwell on the discovery of the X when both 
of the pair-produced Xs decay to Z+j. For the purposes 
of this paper, we will restrict ourselves to looking at lep- 
tonic decay modes of the Z as it has considerably less 
background, but a full collider study can combine vari- 
ous channels and extend the reach in mx- 

We will then use the value of mx extracted from this 
analysis in order to discover the Higgs in the events where 
one of the Xs decays to Z + j and the other to h + j. We 
will focus on a scenario with a moderately heavy Higgs, 
with a dominant decay mode to W~^W~, although dis- 
covery of a light Higgs through X production should be 
competitive with Higgs production from the SM as well 
(possibly utilizing the recent search methods involving 
boosted final states [3l|, [S^ ) . Once again, we will limit 
ourselves to leptonic decays of the Z as well as the Ws, 
but in a more detailed analysis, several channels can be 
combined to extend the discovery reach. 

Let us start with the Tevatron bounds on mx- The 
strongest constraint arises from a recent analysis of WZ 
production [sl], searching for events with exactly three 
final state leptons. In our benchmark model, events with 
pair-produced X decaying to Z, h and jets, followed by 
h WW can be picked up by this search, as well as 
ZZ plus jet final states when both Z decay leptonically 



but one lepton fails to be identified. We implemented 
the cuts described in Ref. [33j . which include requir- 
ing the presence of an e+e~ or /i^/i^ pair with invari- 
ant mass in the interval [86 GeV, 106 GeV] as well as 
PT,£i > 20 GeV, PT.{i2,ts) > 10 GeV and > 25 GeV. 
We find that for mx > 300 GcV, the number of events 
from X production after the selection cuts falls within 
one standard deviation of the SM expectation. The WW 
final state poses a weaker constraint, because in this fi- 
nal state a veto on hard jets is imposed in order to re- 
duce the ti background Q- For "m-x > 300 GeV the 
X-production cross-section is only a fraction of the ti 
cross-section, so there is no additional constraint from 
searches for tt either. Finally, for the same range in mx, 
the Higgs production from X decays is significantly below 
the SM Higgs production cross-section for the same m^, 
so Higgs searches also do not give additional constraints. 
In order to stay consistent with these constraints, we 
will concentrate on the rest of this section on two mass 
points, mx = 300 GcV and mx ~ 550 GeV, both with 
mh = 200 GeV. 

In our analysis, we generate 3 x 10^ parton level events 
(at 14- TeV LHC running) for each value of mx using 
the user-mode of MadGraph [3^ and CTEQ6 parton dis- 
tribution functions [SQl. We decay the X, h, W and Z 
particles using BRIDGE [36[ in order for the angular distri- 
butions of the final-state leptons and jets to be accurate. 
In the signal sample we allow all possible decays of X, 
h, W and Z in order to take full account of the com- 
binatoric background, while in the background samples 
we force leptonic decays (including t's) since hadroni- 
cally decaying background events will not pass our se- 
lection cuts. The hadronization and detector simulation 
arc done with the PYTHIA [sj and PGS [2l[ interface in 
MadGraph, with the default CMS parameter set. This pa- 
rameter set uses a tracker and muon system rj coverage 
up to 2.4, and a minimum lepton px of 5 GeV. Cone jets 
with Ai? = 0.5 are used, note however that our anal- 
ysis is inclusive and we do not place any cuts on jets. 
The PGS default algorithms arc used for lepton isolation 
as well as other details of detector simulation. As the 
background for the X search we generated a matched 
ZZ + jets sample containing 53494 events after match- 
ing, using the MLM matching in MadGraph, and the same 
tools as were used for the signal. For the Higgs search we 
generated a matched tiZ -I- jet sample with 36120 events 
after matching in the same way. 



1. Discovering the X. 

In order to discover the X, we focus on the ZZ final 
state, where both Zs decay leptonically. For SM back- 
grounds, we have generated an MLM-matched sample 
with ZZ, ZZ + j, and ZZ + 2j. As our event selec- 
tion criteria, we demand that an event contains two (dis- 
tinct) Z candidates, where a Z candidate is defined as 
an e'^e~ or a pair with an invariant mass within 
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FIG. 11: Distribution of the average Z+]et pair mass for 
signal {mx ~ 300 GeV) and background in 10 fb^^ of LHC 
data (14 TeV). 



FIG. 12: Distribution of the average Z+jet pair mass for 
signal {mx ~ 550 GeV) and background in 10 fb^^ of LHC 
data (14 TeV). 



5 GcV of mz- We then pair the two Z candidates with 
the two hardest jets in the event (which in the case of 
signal are expected to come from the partons of the X 
decays) and retain the pairing where the Zj pair masses 
are closer to each other. We then plot the average pair 
mass mzj, where the peak from the X is clearly visible 
and distinguishable from the ZZ + jets background. For 
nix = 300 GeV, where the branching fraction X Zj is 
0.76, the cross-section before (after) selection is 36.8 pb 
(50.5 fb). Similarly, before (after) selection we obtained 
1.43 pb (1.32 fb) for mx = 550 GeV, where the branching 
fraction X — >■ Zj is 0.57. Lcptonic branching fractions 
and selection cuts reduce the SM background cross sec- 
tion from 11.4 pb down to 5.1 fb. The results for the two 
mass points and for 10 fb~^ are plotted in Fig. [TT] and 
Fig. [T^ Since the background peaks towards low values 
of mzj, a cut on mzj can improve signal significance. 
For the mx = 550 GeV case we use mzj > 430 GeV as 
a selection cut. While a similar cut can also be used for 
the mx = 300 GeV case, signal is already much larger 
than background in this case and therefore a cut on mzj 
is not essential. With this additional cut, the signal cross 
section becomes 0.85 fb while background is reduced to 
0.074 fb. With an integrated luminosity of 10 fb~^, this 
translates to an average 9.2 events with a background 
expectation of 0.74 events. Using Poisson statistics, this 
corresponds to a probability of 9.4 x 10~®, equivalent to 
more than a 5a upward fluctuation in a Gaussian dis- 
tribution. We conclude that mx up to « 550 GeV is 
discoverable with 10 fb""^ at = 14 TeV. 



2. Discovering the Higgs 

In order to discover the Higgs, we focus on the Zh final 
state, where the Z as well as the Ws from the Higgs decay 
to leptons. More concretely, our event selection criteria 



are: 

• The event contains two positively and two nega- 
tively charged leptons, with exactly one Z candi- 
date (as defined above). 

• The Z candidate and (at least) one of the two hard- 
est jets has an invariant mass within 90 GeV of mx 
as determined in the previous section. 

For mx ~ 300 GeV the selection reduced the signal cross- 
section from 36.8 pb to 22.2 fb, for mx = 550 GeV from 
1.43 pb to 0.87 fb. With the first selection criterion, the 
dominant background is tiZ -I- jet(s). At parton level, 
we generated an MLM-matched sample with up to one 
extra parton, i.e., tiZ and tiZ + j. The second selection 
criterion, which uses the value for mx obtained with the 
search strategy described in the previous subsection, then 
further reduces the background such that we are essen- 
tially left with a pure signal sample. Leptonic branch- 
ing fractions and the selection cuts reduce the back- 
ground cross-section from 0.61 pb to 0.40 fb (0.083 fb), 
applying the two selection criteria for mx = 300 GeV 
{mx ~ 550 GeV). For both mass points, this corre- 
sponds to discovery level statistical significance. Using 
Poisson statistics in the heavy mass case, the probabil- 
ity for a background fluctuation to mimic the signal is 
2.4 X 10~^, equivalent to more than 5(t in a gaussian dis- 
tribution. 

We then identify the two leptons which do not belong 
to the Z candidate, and form the transverse mass variable 
Mt,ww as follows: 

^T,WW = (EtA+I- +E'T^yyf - {pt,1+1- +Fr,miss)^,(16) 

where 

Et,^u = PT,miss + mf+i- . (17) 
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rrih ~ 200 GeV) and background with a luminosity of 10 fb^^. 
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unification. We have identified the characteristic features 
of such models. New colored particles at the TeV scale 
are ubiquitous, which can be produced at the LHC. The 
colored particles allow couplings such that they domi- 
nantly decay to a Higgs and a jet, within or outside the 
LHC detector depending on the size of the couplings. 
The former possibility gives rise to a new channel for 
Higgs production, while the latter leads to spectacular R- 
hadron signals. In order to study these interesting char- 
acteristic collider signatures of WIMP DM and unifica- 
tion, we have chosen the model with the simplest matter 
content as a benchmark. The benchmark model contains 
two generations of an SU(2)L-doublet color-triplet parti- 
cle X that decay via Yukawa-type couplings to the SM. In 
particular, the final states contain Ws, Zs and Higgses, 
as well as down- type quarks. 

We have then investigated the constraints from fla- 
vor bounds on the size of the Yukawa-type coupling that 
leads to the X decay. We showed that there is a range 
where the X can decay promptly, or can be long-lived 
(stable on collider time scales). We have explored each 
of these possibilities in turn, showing that in the case 
of the long-lived X, R-hadrons can be discovered at the 
early LHC (7 TeV) up to mx = 650 GeV, and past 1 TeV 
with a 14- TeV running. 

In the case where the X particles decay promptly, 
a large number of Higgs bosons are produced through 
the decays, which, depending on the X mass, can be 
discovered with less luminosity than would be possible 
from SM Higgs production. We have shown that, for 
mx < 550 GeV and with 10 fb"^ of data at the LHC, 
we can discover the X itself in the leptonic ZZ final 
state as well as the Higgs bosons from X decays in the 
leptonic WW final state for a benchmark Higgs mass 
mh = 200 GeV. 



Note that Et^up is only an approximation to the true 
transverse energy of the neutrino system 

-^T,i/P,tmc — Pt.vv + "^^P- (18) 

Although this approximation only becomes exact when 
the W are produced at threshold, for a 200-GeV Higgs 
the M'Y^ww distribution still peaks near the Higgs mass. 
We plot the results for the mx = 300 GeV in Fig. [T^ and 
for mx = 550 GeV in Fig. [HI 

VI. CONCLUDING REMARKS 

The minimal way to incorporate a WIMP DM candi- 
date is as the neutral component of an electroweak triplet 
with zero hypcrchargc. We have looked for possible ex- 
tensions of the SM that contain such a triplet as well as 
additional matter fields as necessary for gauge coupling 



While in this paper our goal was to focus only on sig- 
nals that are clean and have little background, these stud- 
ies can be significantly expanded in a more dedicated 
collider search. In particular, semi-leptonic decay chan- 
nels can be combined with the fully leptonic ones to in- 
crease the reach. The discovery potential of a light Higgs 
should be enhanced as well, especially utilizing the re- 
cently introduced search methods relying on boosted final 
states. Finally, models other than the benchmark model 
we have chosen can be studied for qualitatively different 
final states. For example, while final states with up-type 
quarks (in particular the top quark) arc rare in the bench- 
mark model, other models can give rise to a large number 
of tops produced from the decays of the new physics (in 
addition to the two signals we discussed in this paper). 
Another interesting problem is to study the case where 
X decays within the detector with a displaced vertex or 
with a macroscopic length. These questions will be fur- 
ther explored in future work. 
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Appendix A: Forbidding proton decay 

In Sec. ini we did not require the unification scale to 
be high enough to suppress proton decay. For example, 
in the benchmark V + 2X model, the unification scale is 
Mgut ^ 10^^ GeV, significantly lower than the standard 
GUT scale 10^^ GeV. Therefore, it is essential to 
discuss how proton decay can be avoided in principle, 
even though it has no phenomenological relevance at the 
TeV scale. 

A robust way to avoid the proton decay problem is 
simply to forbid it by a symmetry. For example, we can 
consider "baryon triality" [s^, (p — >■ e'^^^*"^^* V, where 
and arc the baryon number and hypercharge of 
a field 0, respectively. Baryon triality implies that the 
baryon number can only be violated in units of 3. Thus, 
it will absolutely forbid not only proton decay but also 
neutron-antineutron oscillation [39|. Alternatively, for 
forbidden proton decay only, one could impose "quark 
parity" , q,d'^,u'^ — >■ —q, —d'^, —u'^. 

Clearly, such symmetries necessarily treat quarks and 
leptons differently, so the simplest possibility of embed- 
ding quarks and leptons into unified GUT multiplets 
(e.g. 5 10 for SU(5)) does not work. Instead, quarks 
and leptons must come from separate multiplets, e.g.: 

where and 5^ carry the same baryon triality (or quark 
parity) quantum numbers as d'^ and £, respectively. Be- 
fore discussing how to promote such incomplete multi- 
plets to fully GUT-invariant multiplets, we would like to 
point out that splitting SM matter fields from their heavy 
"GUT partners" is plausible in the sense that there is al- 
ready a split multiplet, i.e. the Higgs doublet, and what- 
ever mechanism that splits the Higgs from its GUT part- 
ner could also work for the SM matter fields. Moreover, 
separating quarks and leptons allows us to trivially in- 
corporate non-unified mass relations for the 1st and 2nd 
generations (i.e. me/md, m^/mc, ^ 1 at the unification 
scale). 

Now, let us discuss how the above "incomplete" multi- 
plets can be compatible with a GUT symmetry. A simple 
and plausible way to do so is to copy the mechanism [40[ 
of "incomplete" multiplets nature already has in the 



low-energy QCD. QCD undergoes a symmetry breaking 
from G = SU(3)l ® SU(3)r down to H = SU(3)l©r. 
Here, the low-mass hadrons form "incomplete" multi- 
plets" , i.e. multiplets of the unbroken subgroup H, not 
those of the full group G. Of course, the theory must be 
invariant under G, which is accomplished by nonlinear 
realization in terms of the Nambu-Goldstonc boson 
field S transforming as S — > gT,h~^ with g ^ G and 
h £ H. This permits us to promote any multiplet of H, 
(pH, to a full multiplet of G, 4>h, as (pa = Ticpn- 

By analogy, let us consider the following GUT scenario. 
Imagine a new confining strong dynamics with a "fiavor 
symmetry" G = SU(5) which undergoes "chiral symme- 
try breaking" G ^ H with H = SU(3) ® SU(2) ® U(l). 
G is also weakly gauged, providing the SM gauge group 
with a single unified coupling at the G confinement scale 
(which therefore can be called the GUT scale). We then 
imagine that the split quark and Icpton multiplets as well 
as the Higgs doublet are "hadrons" of the new strong 
dynamics which form multiplets of iJ, with appropri- 
ate quark parity or baryon triality, etc. They can be 
promoted to full G multiplets by nonlinear realization. 
(If we further integrate in the "p mesons" into this pic- 
ture by employing "hidden local symmetry" (4^ , then we 
essentially obtain the 2-site moose model considered in 
Ref. H.) 

Constructing an explicit, UV-complete 4D gauge the- 
ory that realizes this scenario is beyond the scope of this 
paper, but it is straightforward to construct a 5D real- 
ization via the AdS / GET correspondence [3l , as is done 
in Ref. [l^. The sim ples t setup would be the Randall- 
Sundrum framework [46| , i.e., where we take the UV 
boundary and the bulk to be G-symmetric while the IR 
boundary to be only i?-symmetric, by putting a G gauge 
field in the bulk with the boundary condition that the 
G/ H gauge bosons vanish at the IR boundary. The split 
matter and Higgs multiplets can then be put either in the 
bulk with the appropriate boundary conditions to project 
out the unwanted "G-partner" states, or simply at the IR 
boundary as multiplets of H . With an appropriate sym- 
metry such as quark parity or baryon triality, this setup 
solves the proton decay problem, but since the scale as- 
sociated with the IR boundary is the GUT scale, there 
are no observable consequences of the Kaluza-Klein exci- 
tations. (There are, however, RS GUT models with the 
TeV-scale IR brane, with split multiplets and a discrete 
symmetry to forbid proton decay, which do have exper- 
imental consequences at the TeV scale [13]. The main 
difference between those models and ours is that their 
physics above the TeV scale is a strong conformal dy- 
namics (or a 5D theory via AdS/CET) while we assume 
a perturbative 4D physics up to the GUT scale.) 



Appendix B: Effects of higher-dimensional operators 

Here we would like to show that adding nonrenormal- 
izable interactions to our Lagrangian ([1} does not alter 
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our results. In particular, one might worry that a higher- 
dimensional operator might lead to a prompt decay of the 
X even when A = 0, destroying the R-hadron scenario. 
Fortunately, this is not the case. The leading nonrenor- 
malizable interaction that can let X decay is 



decay length due to this operator is at least 



mxv 



167r(10" GeV)2 



3 m ■ 



1 TeV 
mx 



(B2) 



j{Xm*){qH*), 



(Bl) 



where A is some high scale. By our assumption, there 
is no new threshold between the TeV scale and the uni- 
fication scale, so A must be at least ~ 10"'^^ GeV. To 
avoid a suppression from three-body phase space, one of 
the Higgs fields can be put to its VEV. Therefore, the X 



Thus even the most conservative estimate gives a decay 
length comparable to the dimensions of the LHC detec- 
tors. Therefore, higher-dimensional operators do not up- 
set our conclusions. Since this X decay is still prompt on 
a cosmological time scale, setting A = in the coupling 
(jl]) is actually allowed cosmologically, as X can decay via 
the above operator. 
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